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Abstract. Time evolution of the expectation values of various dynamical operators of the harmonic 
oscillator with dissipation is analitically obtained within the framework of the Lindblad theory for 
open quantum systems. We deduce the density matrix of the damped harmonic oscillator from 
the solution of the Fokker-Planck equation for the coherent state representation, obtained from the 
master equation for the density operator. The Fokker-Planck equation for the Wigner distribution 
function, subject to either the Gaussian type or the <5-function type of initial conditions, is also 
solved by using the Wang-Uhlenbeck method. The obtained Wigner functions are two-dimensional 
Gaussians with different widths. 



1 Introduction 



In the last two decades, the problem of dissipation in quantum mechanics, i.e. the consistent 
description of open quantum systems, was investigated by various authors [1-6]. Because 
dissipative processes imply irreversibility and, therefore, a preferred direction in time, it 
is generally thought that quantum dynamical semigroups are the basic tools to introduce 
dissipation in quantum mechanics. The most general form of the generators of such semi- 
groups was given by Lindblad [7-9]. This formalism has been studied for the case of damped 
harmonic oscillators [8, 10, 11] and applied to various physical phenomena, for instance, the 
damping of collective modes in deep inelastic collisions in nuclear physics [12, 13] and the 
interaction of a two- level atom with the electromagnetic field [14]. Recently [15], a family 
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of master equations, constructed in the form of Lindblad generators, was proposed for local 
ohmic quantum dissipation. 

This paper, dealing with the damping of the harmonic oscillator within the Lindblad 
theory for open quantum systems, is concerned with the time evolution of various dynamical 
operators involved in the master and Fokker-Planck equations, in particular with the time 
development of the density matrix. In [16] the Lindblad master equation was transformed 
into Fokker-Planck equations for quasiprobability distributions and a comparative study was 
made for the Glauber P, antinormal ordering Q and Wigner W representations. In [17] the 
density matrix of the damped harmonic oscillator was represented by a generating function. 
We shall explore the physical aspects of the Fokker-Planck equation which is the c-number 
equivalent equation to the master equation for the density operator. Generally the master 
equation gains considerably in clarity if it is represented in terms of the Wigner distribution 
function which satisfies the Fokker-Planck equation. It is worth mentioning that these master 
and Fokker-Planck equations agree in form with the corresponding equations formulated in 
quantum optics [18-25]. 

The content of the paper is arranged as follows. In Sec. 2 we review the derivation of the 
master equation of the harmonic oscillator. In order to get an insight into physical meanings 
of this equation, we first split it up into several equations satisfied by the expectation values 
of dynamical operators involved in the master equation. These equations are then solved 
analytically. In Sec. 3 we transform the master equation into the Fokker-Planck equation by 
means of the well-known methods [4,26-30]. We extract the density matrix with the help of 
the solution of the Fokker-Planck equation for the coherent state representation. Then the 
Fokker-Planck equation for the Wigner distribution, subject to either the Gaussian type or 
the ^-function type of initial conditions, is solved by the Wang-Uhlenbeck method. Finally, 
conclusions are given in Sec. 4. 



2 Master equation for the damped harmonic oscillator 

The rigorous formulation for introducing the dissipation into a quantum mechanical system 
is that of quantum dynamical semigroups [2, 3, 7]. According to the axiomatic theory of 
Lindblad [7, 9] , the usual von Neumann-Liouville equation ruling the time evolution of closed 
quantum systems is replaced in the case of open systems by the following equation for the 
density operator p: 

^ = i(*.(p». (2.1) 

Here, $ t denotes the dynamical semigroup describing the irreversible time evolution of the 
open system in the Schrodinger representation and L the infinitesimal generator of the 
dynamical semigroup $ t . Using the structural theorem of Lindblad [7] which gives the most 
general form of the bounded, completely dissipative Liouville operator L, we obtain the 
explicit form of the most general time- homogeneous quantum mechanical Markovian master 
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equation: 



= L{p{t)) = -UH,p(t)} + ^£([^),^] + [V^p{t)V]\). (2.2) 



dp(t) 

dt Kr ^" ^"'^"/J ■ ^ 
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Here H is the Hamiltonian of the system. The operators Vj and Vj are bounded operators 
on the Hilbert space of the Hamiltonian. 

We should like to mention that the Markovian master equations found in the literature 
are of this form after some rearrangement of terms, even for unbounded Liouville operators. 
In this connection we assume that the general form of the master equation given by (2.2) is 
also valid for unbounded Liouville operators. 

In this paper we impose a simple condition to the operators H,Vj,Vj that they are 
functions of the basic observables q and p of the one-dimensional quantum mechanical system 
of such kind that the obtained model is exactly solvable. This condition implies [8] that Vj 
are at most first degree polynomials in p and q and H is at most a second degree polynomial 
in p and q. Because in the linear space of the first degree polynomials in p and q the operators 
p and q give a basis, there exist only two C-linear independent operators Vi, V 2 which can 
be written in the form 

V j = a jP + b j q, j = 1,2, (2.3) 

with dj,bj complex numbers [8]. The constant term is omitted because its contribution to 
the generator L is equivalent to terms in H linear in p and q which for simplicity are assumed 
to be zero. Then the harmonic oscillator Hamiltonian H is chosen of the form 

H = H + ^(pq + qp), H = ±-f + ^fq\ (2.4) 

With these choices and introducing the annihilation and creation operators via the relations 



we have H = hu(a^a +1/2) and the Markovian master equation takes the form 

^j- = -(D 1 + fi)(a^a^ p — a) pa)) H — {D x — p){pa)a) — a) pa)) 
dt 2 2 

+ -(D 2 + A + iu)(apa^ — a^ap) H — (D 2 — A — iu)(a^pa — paa^) + H.c, (2.6) 
2 2 



where 



D 1 = UmcuD qq -^ + 2iD pq ), D 2 = \{mwD m + ^). (2.7) 



Here we used the notations: 



D m = 9 J2 M » D pp = 9 J2 N » D pq = D qp = ~9 Re J2 a j b J> X = ~ Im J2 a j b J> 

Z 3=1,2 Z j=l,2 Z j=l,2 3=1,2 

(2.8) 
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where D^,, D qq and D pq are the diffusion coefficients and A the friction constant. They satisfy 
the following fundamental constraints [10]: 

i) D pp > 0, ii) D qq > 0, iii) D pp D qq - D 2 q > AV/4. (2.9) 

In the particular case when the asymptotic state is a Gibbs state 

p G (oo) = e-^/Tre-^, (2.10) 

these coefficients reduce to 

A + /i hu X — fj, h hcu 

Am = nmuj coth——, D aa = coth — — , D va — 0, (2.11) 

pp 2 2kT' qq 2 muj 2kT' K J 

where T is the temperature of the thermal bath. 

In the literature, master equations of the type (2.6) are encountered in concrete theoretical 
models for the description of different physical phenomena in quantum optics [18-24], in 
treatments of the damping of collective modes in deep inelastic collisions of heavy ions [31- 
34] or in the quantum mechanical description of the dissipation for the one-dimensional 
harmonic oscillator [4, 6, 28, 29]. A classification of these equations, whether they satisfy or 
not the fundamental constraints (2.9), was given in [17]. 

The meaning of the master equation becomes clear when we transform it into equations 
satisfied by various expectation values of operators involved in the master equation, < A >= 
Tr[p(t)A], where A is an operator composed of the creation and annihilation operators. 
Multiplying both sides of (2.6) by a and taking throughout the trace, we get 

-f- < a >= -(A + iuj) < a > +/i < a ] > . (2.12) 
at 



Similarly, the equation for < a) > is given by 



Y < a) >= -(A - iuj) < a ] > +/i < a > . (2.13) 

In the absence of the second term on the right-hand side of (2.12) and (2.13), the two equa- 
tions represent independently a simple equation of oscillation with damping. By coupling 
(2.12) to (2.13) we get a second order differential equation for either < a > or < a) >. For 
example, we obtain: 

d 2 d 

— < a > +2A— < a > +(A 2 + cu 2 - /i 2 ) < a >= 0, (2.14) 

LLL (JjL 

which is the equation of motion for Brownian motion of a classical oscillator, but without the 
term corresponding to random process. Because of the vanishing terms on the right-hand 
side, we may equally state that (2.14) is the equation of motion with zero expectation value 
of the random process. In the study of the Brownian motion of a classical oscillator, one 
replaces the second-order differential equation of motion of the type (2.14) by two equivalent 
first-order differential equations [35] which are precisely the Langevin equations. Accordingly, 
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we may say that (2.12) and (2.13) are the Langevin equations corresponding to (2.14), but 
without the random process term. The integration of (2.14) is straightforward. There are 
two cases: a)// > uj (overdamped) and b)/x < uj (underdamped). In the case a) with the 
notation v 2 = ji 2 — uj 2 , we obtain: 

< a {t) >= e~ xt [< a(0) > (coshz/t - i—wahut) + - < a f (0) > sinhi/i]. (2.15a) 
In the case b) with the notation Vl 2 = cj 2 — /i 2 , we obtain: 

< a (t) >= e ~ A '[< a (0) > (cos fit — i— sin nt) + ^ < a f (0) > sin fit]. (2.15b) 

The expression for < a*(£) > can be obtained simply by taking the complex conjugate of the 
right-hand side of (2.15). 

For the computation of quantal fluctuations of the coordinate and momentum of the 
harmonic oscillator, we need the expectation values of quadratic operators, such as a +2 ,a 2 
or a)a. The dynamical behaviour of these operators can be well surveyed by deriving the 
equations satisfied by their expectation values. By following the same procedure as before, 
employed in the derivation of (2.12), we find: 

4- < a 2 > +2(A + iu) < a 2 >= 1\i < a ] a > +D X + fi, (2.16) 
at 

4 < a +2 > +2(A - iu) < a +2 >= 2/x < a ] a > +D\ + fi, (2.17) 
at 

4- < a ] a > +2A < a f a >= /i(< a +2 > + < a 2 >) + D 2 - A. (2.18) 
at 

The solutions of these equations are readily obtained by transforming them into two differ- 
ential equations satisfied by the sum and the difference of two quadratic operators, < a 2 > 
and < a +2 >. We obtain: 

\^{< a 2 > + < a +2 >) + A^(< a 2 > + < a +2 >) + (A 2 + uj 2 - /i 2 )(< a 2 > + < a +2 >) 
Adt 2 dt 

= I[(A + pi)mu;D qq - (A - /i)^ + 2uD m \ = D, (2.19) 
n muj 

~(< a 2 > - < a +2 >) + A(< a 2 > - < a +2 >) + iu(< a 2 > + < a +2 >) = 2i^. (2.20) 

jL (XL ft 

The solution of (2.19) is straightforward and with the help of which both equations (2.18) 
and (2.20) can be immediately solved. We find: 

< a 2 >= e- 2 ^[(l - £)Cie™ + (1 + i-)C 2 e~™ - i-C 3 ] + ££f % + (2.21a) 
v v uj 2\{\ l — v l ) TiX 

< a^a >= e'^C^ - C 2 e~ 2vt ) + C 3 ] + ^{^~ 2 + D 2 - A), (2.22a) 

V LA A — V z 



5 



< a)a >- 



for the overdamped case p > u and 

< a 2 >= e-^KC! + *C 2 £) cos2^ + (C 2 - ic£) sin2fii - ^C 3 ] + ^"'ol + & 

a; 2A(A Z + ii J ftA 

(2.21b) 

" 2At [^(Ci sin2fit - C 2 cos2fit) + C 3 ] + ^(^^ + D 2 - A), (2.22b) 

for the underdamped case uj > p. The expression for < a +2 > can be obtained by taking 
the complex conjugates of (2.21a) and (2.21b). Here, Ci, 62,63 are the integral constants 
depending on the initial expectation values of the operators under consideration. In partic- 
ular, if D qq = D pq = and p = A we obtain for the underdamped case uj > p the equations 
written by Jang [31] for the model on nuclear dynamics based on the second RPA at finite 
temperature. For time t — > 00, we see from (2.22) that 

<B ' B>= ^ + S'-^ + '- A »- (2 ' 23) 

In the particular case when the asymptotic state is a Gibbs state (2.10), we get 

< ata >= -(coth — - 1) = (exp — - l)- 1 =< n >, (2.24) 

which is the Bose distribution. This means that the expectation value of the number operator 
goes to the average thermal-phonon number at infinity of time. From the identity 

00 

< a^a >= ^2 m < m\p(t)\m > (2.25) 

m=0 

it follows 

< n > m 

< m\p(po)\m >= — -. (2.26) 

irv i\ (l+<n>) m+1 v 1 

In deriving this formula, we have made use of the identity J2m=o mxm = x /(l ~ X Y '■ The 
expression (2.26) shows that in the considered particular case the density matrix reaches its 
thermal equilibrium - the Bose-Einstein distribution, whatever the initial distribution of the 
density matrix may be. When the initial density matrix < m\p(0)\m > is represented by a 
distribution of the form N m /(1 + iV) m+1 , where N stands for the average phonon number, 
the relation (2.25) implies that < a^a >= N. When the initial density matrix is characterized 
by a distribution of the form 

— t N m e~ N , (2.27) 
ml 

(2.25) implies that < a^a > becomes also N. Eq. (2.27) is nothing but a Poisson distribution. 
If the initial density matrix is represented by a Kronecker delta 5 ms , we see from (2.25) that 

< a)a >= s, which corresponds to the initial s-phonon state. 

The physical observables of the harmonic oscillator can be obtained from the expectation 
values of polynomials of the annihilation and creation operators. So, for the position and 
momentum operators q and p via the relations (2.5), we can evaluate either the second 
moments or variances (fluctuations), by making use of the results (2.15), (2.21), (2.22). 
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3 Fokker-Planck equations 



One useful way to study the consequences of the master equation (2.6) for the density oper- 
ator of the one-dimensional damped harmonic oscillator is to transform it into more familiar 
forms, such as the equations for the c-number quasiprobability distributions Glauber P, anti- 
normal ordering Q and Wigner W associated with the density operator [16]. In this case the 
resulting differential equations of the Fokker-Planck type for the distribution functions can 
be solved by standard methods [26-29] employed in quantum optics and observables directly 
calculated as correlations of these distribution functions. 



3.1 Calculation of the density matrix from the Fokker-Planck equa- 
tion 

The Fokker-Planck equation, obtained from the master equation and satisfied by the Wigner 
distribution function W(a,a*,t), where a is a complex variable, has the form [16]: 

dW(a a* t) d d 

' = -{ — [-(X + ioo)a + fia*} + — — [-(A - iuo)a* + fia]}W{a, a*, t) 

at da da* 

When we substitute the P representation function P(a,a*,t) for W(a,a*,t) and the coeffi- 
cients D\ +/i for Di, D 2 — A for D 2 in the above equation, we get the Fokker-Planck equation 
for the coherent representation [16]. 

The Fokker-Planck equation for the P representation, subject to the initial condition 

P(a,a*,0) =5(a-a )5(a* -a* ), (3.2) 

where ao is the initial value of a can be solved [30, 36, 37] and the solution (Green function) 
for the P representation is found to be 

2 1 

P(a,a*,t) = — exp{- [a 22 (a - a ) 2 + a u (a* - a^) 2 - 2a 12 \a - a \ 2 } , (3.3) 

7i7det<7(f) det(7 W 



where 

<T %j( t )= Yj [Sisdjr-b is (t)b jr (t)]a sr (oo). (3.4) 

s,r=l,2 

The function d Q and its complex conjugate, which are still functions of time, are given by 

«o = b n (t)a + b 12 (t)a* . (3.5) 
The functions by obey the equations 

ks = Ci 3 h 3s (3-6) 
i=i,2 
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with the initial conditions bj S (0) = 5j S and cr(oo) is determined by 

CV(oo) + a(oo)C T = Q p , (3.7) 

where 

C= { -p A-iJ' ^ = {d 2 -X Dl + n)- (3 ' 8) 

We get 

b n =b* 22 = e~ xt (cos ttt-i^ sin tit), b 12 = b 21 = sin fit, (3.9) 

with fi 2 = cj 2 — yU 2 . While the functions au,a 22 and a are complex with <7n = a 22 , the 
functions deter (t) and o"i 2 are real. 

The solution of the Fokker-Planck equation has been written down providing the diffusion 
matrix Q p is positive definite. However, the diffusion matrix in the Glauber P representation 
is not, in general, positive definite. If the P distribution does not exist as a well-behaved 
function, the so-called generalized P distributions can be taken that are well-behaved, normal 
ordering functions 



In the coherent representation [39, 40] the density operator p(t) is expressed by 

p[t) = J P{a,a*,t)\a><a\d 2 a, (3.10) 

where d 2 a = d{Rea)d(lma) and \a > is the coherent state. The matrix element of p(t) in 
the n quantum number representation is obtained by multiplying (3.10) on the left by < m\ 
and on the right by \n > . By making use of the well-known relation 

1 00 a n 

\a >= exp(--|a| ) ^ -^=|n >, (3.11) 
I n =o v n\ 

we get 

<m\p(t)\n>= J: / a n a* m P(a,a*,t)exp(-|a| 2 )rf 2 a. (3.12) 

ym\n\ J 

Upon introducing the explicit form (3.3) for P(a,a*,t) into (3.12), we obtain the desired 
density matrix for the initial coherent state. However, due to the powers of complex variables 
a and a* in the integrand, the practical evaluation of the integral in (3.12) is not an easy 
task. Instead, we use the method of generating function [31] which allows us to transform 
(3.12) into a multiple-differential form. When we define a generating function F(x,y,t) by 
the integral 

F(x,y,t) — J P(a, a* ,t) exp(— |a| 2 + xa + ya*)d 2 a, (3.13) 
we see that the density matrix is related to the generating function by 

< m\p(t)\n >= -yl=(^-r(^-rF(x,y,t)\ x=y=0 . (3.14) 
\Jm\n\ ox dy 
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Since the P representation is in a Gaussian form, the right-hand side of (3.13) can be 
evaluated analitically by making use of the identity 

f ( i i2 , 7 i * , 2 . t *2x,2 71 abc + b 2 f + c 2 e 

exp -a z + bz + cz + ez + fz )d z = = = exp — , (3.15) 

J V a — 4ej cr — 4ej 

which is convergent for Re a > \e* + f\, while b, c may be arbitrary. We find: 
F = -j=exp{xy - ^[a n (x - a* ) 2 + a 22 (y - a ) 2 - 2(a 12 + 2)(x - a* )(y - a )]}, (3.16) 



where A = d — 4(cri 2 + 1), d = deter = aua 22 — o\ 2 - A formula for the density matrix can 
be written down by applying the relation (3.14) to the generating function (3.16). We get 

< m\p(t)\n >= 2 r =^ exp{[--[a 22 al + a^af - 2(a 12 + 2)\a \ 2 }} 
J\A\ A 

( — 1 \rii+n 2 r)m+n-2(n 1 +n2+n- s ) 

X E , , \, 7r=- TT7 z rrE, (3.17) 

m,nt^3=o «i!n 2 !n 3 !(m - 2m - n 3 )!(n - 2n 2 - n 3 )\ 

where 

<7>22 2 (d - 2a 12 ) n 3[a n aS - (m 2 + 2)a p- 2ni - n3 [<x 22 a - 2(<r 12 + 2)al\ n ~ 2n ^ 



E = 



j^m+n-(n 1 +n2+n :i ) 



The expression (3.17) is the density matrix corresponding to the initial coherent state. At 
time t — 0, the functions mi, cr 22 and <j\ 2 vanish and a goes to a . In this case the density 
matrix reduces to 

< m\p(0)\n >= -jl=a^a* m ex P (-\a \ 2 ), (3.18) 
ym\n\ 

which is the initial Glauber packet. For the diagonal case the initial density matrix becomes 
the Poisson distribution. At infinity of time, the density matrix (3.17) goes to the Bose- 
Einstein distribution 

< n > m 

< m\p(oo)\n >= (1+<n>)m+1 <W (3.19) 



3.2 Wigner distribution function 

The Fokker-Planck equation (3.1) can also be written in terms of real coordinates x\ and x 2 
(or the averaged position and momentum coordinates of the harmonic oscillator) defined by 
a = xi + ix 2 , a* = Xi — ix 2 , as follows: 

f - E/^W^ EC^ (3.20) 

where 

\ uj A + /J' * ~h\ D pq Djmuj- [6 - ZL) 
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Since the drift coefficients are linear in the variables x\ and x 2 and the diffusion coefficients 
are constant with respect to x\ and x 2 , (3.20) describes an Ornstein-Uhlenbeck process 
[35, 41]. Following the method developed by Wang and Uhlenbeck [35], we shall solve this 
Fokker-Planck equation, subject to either the wave-packet type or the 5-function type of 
initial conditions. By changing the variables x\ and x 2 of (3.20) via the relations 

Z\ = axi + bx 2 , z 2 = cxi + dx 2 , (3.22) 

(3.20) is transformed into the standard partial differential equation [35] expressed as 

- <-E<* ^E/«*^(^.0. (S.2S) 

In deriving this equation we have put a — c* — (// — b — d — 1, v\ — = —A — iVl 

and 

D n = D* 22 = -*-[(// - iti) 2 mD qq + 2(/i - in)D pq + ^], 
nuj m 

D 12 = D 21 = UmLuD qq + ^D pq + ^). (3.24) 
n uj muj 

1) When the Fokker-Planck equation for the coherent state representation is subject to 
the initial condition 5(a — ao)5(a* — a^), then the use of the relation between the Wigner 
distribution function and P representation 

W(a, a*, t) = - f P(/3, /F, t) exp(-2|a - (3\ 2 )d 2 (3 (3.25) 

7T J 

leads to a Gaussian form for the initial Wigner function. If this Wigner function is expressed 
in terms of x±q and x 2 q ~~ the initial values of x\ and x 2 at t = 0, respectively, then we get 
the expression which corresponds to the initial condition of a wave packet: 

W w (x 1 ,x 2 ,0) = ^W w (a,a*,0) = exp(-2|a-a | 2 ) = ^ exp{-2[(xi-a;io) 2 +(x2-a;2o) 2 ]}- 

(3.26) 

Accordingly, we now look for the solution of the Fokker-Planck equation (3.20) subject to 
(3.26). By changing the variables x\ and x 2 into z\ and z 2 , [z\ — z 2 = z), this initial 
condition is seen to be transformed into 

1 2uj 2 

W w (z, z\ 0) = — exp{ — [q(z - z ) 2 + q*(z* - z*) 2 - \z - z \ 2 }}, (3.27) 

where z is the initial value of z and q = //(// + iVl)/2uo 2 . The solution of (3.23) subject to 
the initial condition (3.27) is found to be 

W w (z,z*,t) = ^^exp{-JL[^ - z e^) 2 + 9l (z* - z* e^) 2 - g 3 \z - z e^\ 2 }}, 

(3.28) 
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where 

B w = gm - \gl, gi = g* 2 = g*e 2 ^ + gi( e 2 ^ - l),g 3 = e" 2 * + ^(1 - e~™). (3.29) 
In terms of real variables X\ and x 2 we have: 

W w (x 1: x 2 , t) = . exp{--^—[(f) w (xi - xif + tp w (x 2 - x 2 ) 2 + Xw(xi - x 1 )(x 2 - x 2 )]}, 

irhwJlBJ IB ™ 



(3.30) 

where 

4>w = gia* 2 + g2a 2 ~ 9s, if> w =9i+92- gs, Xw = + g 2 a) - g 3 (a + a*). (3.31) 

The functions X\ and x 2 , which are also oscillating functions, are given by 

x 1 = e~ xt [x 10 (cosQt + ^ sin fit) + ^20^ sin fit], 

x 2 = e _A< [x 20 (cos fit — ^ sin fit) — ^10^ sin fit]. (3.32) 

2) If the Fokker-Planck equation (3.23) is subject to the 5-function type of initial condi- 
tion, the Wigner distribution function is given by 

W(z,z*,t) = ^exp{-hf 2 (z-z e^ + f 1 (z*-z*e^) 2 -2f 3 \z-z e^\ 2 }}, (3.33) 

7ihuj^\B\ B 

where 

B = hf 2 -fl h = n = — is 2 " 1 * ~ 1), /3 = ^(l-e- 2A '). (3.34) 
In terms of real variables X\ and x 2 we have: 

W(xi,x 2 ,t) = ^^exp{-^[(f) d (x 1 -xi) 2 + ^ d (x 2 -x 2 ) 2 + Xd(x 1 -x 1 )(x 2 -x 2 )]}, (3.35) 

ntujy/\B\ B 

where 

<Pd = ha* 2 + f 2 a 2 - 2/3, Vd = fi + h ~ 2/3, Xd = 2[M* + f 2 a - f 3 (a + a*)}- (3.36) 

So, one gets a 2-dimensional Gaussian distribution with the average values x\ and x 2 and 
the variances (f>d,tpd an d Xd- 

When time t — > 00, X\ and x 2 vanish and we obtain the steady state solution: 

W(x u x 2 ) = 1 = exp[-l 2 (O^(oo)^]. (3.37) 

27r v /det(T^(cx)) 2 iJ=l,2 

The stationary covariance matrix a w (oo) can be determined from the algebraic equation 

Aa w (oo) + a w (oo)A T = Q w . (3.38) 
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We obtain: 




(2A(A + /i) + u 2 )QY + uj 2 QY 2 + 2u(X + /x)Q 
4A(A 2 + u 2 - /I 2 ) 



w 

12 




u 2 QW + (2A(A - /i) + io 2 )Q 2 ^ - 2cu(A - n)Q. 
4A(A 2 + cj 2 - ^ 2 ) 



,w 

12 



(3.39) 




cu(A + l i)QY 1 + cu(A - /i)Q 2 ^ + 2(A 2 
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4A(A 2 + u 2 - /i 2 ) 



4 Conclusions 



Recently we assist to a revival of interest in quantum Brownian motion as a paradigm of 
quantum open systems. There are many motivations. The possibility of preparing systems 
in macroscopic quantum states led to the problems of dissipation in tunneling and of loss 
of quantum coherence (decoherence). These problems are intimately related to the issue of 
quantum-to-classical transition. All of them point the necessity of a better understanding 
of open quantum systems and all requires the extension of the model of quantum Brow- 
nian motion. The Lindblad theory provides a selfconsistent treatment of damping as a 
possible extension of quantum mechanics to open systems. In the present paper we have 
studied the one-dimensional harmonic oscillator with dissipation within the framework of 
this theory. We have carried out a calculation of the expectation values of various dynamical 
operators involved in the master equation, especially the first two moments and the density 
matrix. Generally, the time evolution of the density matrix as well as the expectation values 
of dynamical operators are characterized by complex functions with an oscillating element 
exp(±i^u 2 — fj, 2 t) multiplied by the damping factor exp(— At). We deduced the density ma- 
trix from the solution of the Fokker-Planck equation for the coherent state representation, 
obtained from the master equation for the density operator. For a thermal bath, when the 
asymptotic state is a Gibbs state, a Bose-Einstein distribution results as density matrix. 
The density matrix can be used in various physical applications where a Bosonic degree 
of freedom moving in a harmonic oscillator potential is damped. For example, one needs 
to determine nondiagonal transition elements of the density matrix, if an oscillator is per- 
turbed by a weak electromagnetic field in addition to its coupling to a heat bath. From the 
master equation of the damped quantum oscillator we have also derived the corresponding 
Fokker-Planck equation in the Wigner W representation. The obtained equation describes 
an Ornstein-Uhlenbeck process. By using the Wang-Uhlenbeck method we have solved this 
equation for the Wigner function, subject to either the Gaussian type or the ^-function type 
of initial conditions and showed that the Wigner functions are two-dimensional Gaussians 
with different widths. In a forthcoming paper, by using this Wigner function, it will be 
discussed the entropy of the damped harmonic oscillator. 
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